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ABSTRACT 


Presented  are  the  results  of  the  continuing  study  of  the 
statistical  properties  of  the  envelope  functions  of 
narrowband  reverberation.  Included  are  some  initial 
estimations  of  the  intensity  spectrum  of  the  reverberation 
envelope  processes  generated  by  scattering  1.0  msec  cw 
signals  from  the  rou^,  moving  surface  of  a fresh  water 
lake.  In  addition  the  intensity  spectnjm  of  the  time- 
averaged  envelope  is  estimated.  These  experimental 
estimates  are  compared  to  the  corresponding  theoretical 
quantities  derived  in  earlier  research. 

'\ 


EXPERIMENTAL  ESTIMATES  OF  THE  INTENSITY  SPECTRA 
OF  THE  REVERBERATION  ENVELOPE  AND  ITS  TIME  AVERAGE 

by 

T.  D.  Plemons 


TABLE  OF  CONTENTS 

Page 

ABSTRACT  ill 

I . INTRODUCTION  1 

IT.  REVIEW  OF  THEORY  3 

III.  DESCRIPTION  OF  THE  EXPERIMENTAL  DATA  7 

IV.  THE  FLUCTUATING  CCMPONENT  OF  THE  REVERBERATION  ENVELOPE  9 

V.  ESTIMATION  OF  THE  INTENSITY  SPECTRUM  13 

VI.  RESEARCH  PLANNED  FOR  THE  NEXT  QUARTER  17 

19 


REFERENCES 


t 


I . INTRODUCTION 

Two  previous  progress  reports'*^  in  this  series  present  the 
theoretical  results  obtained  in  the  study  of  the  statistical  properties 
of  the  envelope  functions  of  narrowband  reverberation  processes.  In 
this  report  the  first  experimental  results  of  this  study  are  presented 
with  comparisons  with  the  corresponding  theoretical  expressions. 

Before  proceeding  with  the  discussion  of  the  experimental  data,  a 
review  of  the  theory  will  be  given. 
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II.  REVIEW  OF  THEORY 

The  narrowband  reverberation  process  X(t)  is  of  the  form 

X(t)  = E(t)  cos  1 2nf^t  + i(i(t)j  , ( 

where  4i(t),  and  E(t)  are,  respectively,  the  center  frequency, 
the  phase,  and  envelope  of  X(t).  The  first-  and  second-order 
statistics  of  the  envelope  E(t)  are  the  subject  of  this  study;  in 
particular,  the  covariance  and  spectra  of  E(t)  are  the  desired 
quantities.  Actually  it  is  the  fluctuating  component, 

y(t)  i E(t)  - <E(t)>  , ( 

that  is  of  interest.  Here  ^E(t))  is  the  expected  value  (mean)  of 
E(t).  lYie  covarlemce  of  y(t)  is  (Eq.  (3.9),  Ref.  l) 

K Jti,t2)  H <y(t^)y(t2)> 


fYk/(T)  , T-t^-t^ 


where  f is  the  meeui  intensity  of  the  narrowband  reverberation  X(t) 

and  k (t)  Is  the  envelope  of  the  normalized  covariance  of  X(t). 
o 

Applying  the  Wiener-Khlntchine  theorem  the  dependence  of  the 
associated  intensity  spectrum  of  y(t)  on  the  frequency  f was  derived 
(Eq.  (I4.I),  Ref.  l): 
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with  10  = 2-nf. 


The  next  assumption  was  that  the  covariance  of  the  reverberation 
X(t)  is  equivalent  to  the  autocorrelation  of  the  transmitted  signal, 
S(t).  This  assumption  allowed  us  to  write 

^o(^)  = , (5) 

where  C^(t)  Is  the  envelope  of  the  correlation  function  of  S(t). 
Substituting  this  expression  in  Eq.  (1<)  gave 


W (f)  = 

y 


c/(t) 


cos  (oTdr 


(6) 


Specifying  a transmitted  signal  waveform  determines  C (x).  There 

o 

are  two  cases  of  interest  in  this  study:  frequency  modulated  (JW) 

signals  and  pulsed  cw  signals.  Ibe  calculated  Intensity  spectra  of 
y(t)  corresponding  to  these  signals  are^ 


elsewhere. 


and 


(7) 


TV 
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elsewhere. 


(8) 


Here  T and  are  the  intensities  of  the  narrowband  reverberation 
cw  FM 

processes  generated  by  these  signals,  T the  transmitted  signal  lengths, 
and  W the  bandwidth  of  the  signal.  The  corresponding  covariance 
functions  are 

K (t)  = 5 V (l  - 

y cw  8 cw  \ T / 

= 0 

and 


T < T 


elsewhere , 


(9) 


K - JL  u, 


< T 


elsewhere. 


(10) 


These  results  were  used  to  determine  the  covariance  and  intensity 

spectrum  of  the  envelope  that  is  subjected  to  a time  averaging 

("smoothing"),  which  is  the  case  of  an  incoherent  detector-avera^^er. 

The  fluctuating  component  at  the  output  of  the  detector,  with  an 

averaging  time  T , is 
o 

rt+T  /2 

Y(t)  = — I y(a)da  . (ll) 

o •'t-T  /2 
o 


The  covariance  and  intensity  spectrum  of  Y(t)  are  (Eqs.  (7-2)  and 
(7.3),  Ref.  2) 


and 


o 


WyCf)  = Wy(f) 


sin  ttT  f 
o 

(itT  f)^ 
o 


(12) 


(13) 
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III.  DESCRIPTION  OF  THE  EXPERIMENTAL  DATA 


The  experimental  reverberation  data  were  generated  by  scattering 
sound  from  the  rough,  moving  surface  of  a fresh  water  lake.  The  location 
of  the  experiment  wets  at  the  Lake  Travis  Test  Station  of  Applied  Research 
Laboratories.  The  transmitter  and  receiver  were  located  at  approximately 
the  same  place  8 ft  beneath  the  scatter  siirface  with  their  acoustic  axes 
parallel  and  horizontal;  i.e.,  the  grazing  angle  of  each  acoustic  axis 
was  zero.  A reverberation  ensemble  containing  150  retiims  was  generated 
by  transmitting  cw  pulses  at  a rate  of  1/sec.  The  center  frequency  of 
the  pulse  was  110  kHz  and  the  pulse  duration  was  1 msec. 

The  data  were  recorded  on  analog  magnetic  tape  and  later  converted 
to  digital  samples  to  be  used  for  computer  processing.  Thus  the  quadra- 
ture components  were  sampled  directly^  at  a rate  of  110  kHz.  Figure  1 
shows  computer  plots  of  the  envelopes  of  the  first  three  reverberation 
returns.  These  envelope  functions  are  easily  obtained  by  computing 
the  square  root  of  the  sum  of  the  squares  of  the  quadrature  components. 
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TV.  Tire  FLUCTUATING  COMPONENT  OF  THE  REVEPBEBATTON  ENVELOPE 


The  fluctuating  component  Is  obtained  by  first  computing  the 
ensemble  mean  of  the  envelope  and  then  removing,  through  subtraction, 
this  mean  from  each  envelope  function.  Thus,  we  compute 

y(t)  = E(t)  - <E(t)>  . (ll*) 

The  bar  placed  over  any  svmbol  denotes  the  experimental  estimate  of 
that  quantity.  This  estimate  corresponds  to  Eq.  (2).  The  ensemble 
mean  is  estimated  according  to 


(E(t)) 


150 


= — y 

150  ^ 
1=1 


E^(t) 


(15) 


Thus,  E(t)  is  the  average,  at  a given  time  (range)  t,  of  15O  envelope 
samples . 

In  Fig.  2(c)  the  fluctuating  component  of  the  first  envelope 
function  can  be  seen.  It  was  obtained  by  subtracting  the  envelope 
mean  (average  of  150  returns).  Fig.  2(b),  from  the  envelope  function. 

Fig.  2(a).  It  is  interesting  to  note  that  two  types  of  nonstationarity 
are  observable  in  this  example.  First,  there  is  the  variation,  with 
respect  to  range,  of  the  mean  of  the  envelope  and,  second,  there  is  the 
variation  in  the  Intensity  of  the  fluctuating  component. 

The  fluctuating  component  in  Fig.  2(c)  represents  the  first  of  a 
150-member  ensemble  that  contains  the  second-order  statistics  of  Interest. 
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y,(t)  s E,(t)-<E(t)>  <E(t)> 
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(b)  ENSEMBLE  AVERAGE  OF 
150  REVERBERATION 
ENVELOPE  FUNCTIONS 
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(c)  FLUCTUATING  COMPONENT  IS  THE 
SAMPLE  REVERBERATION  ENVELOPE 
WITH  ENSEMBLE  AVERAGE  REMOVED 


FIGURE  2 

AN  EXAMPLE  OF  THE  FLUCTUATING  COMPONENT 

OF  THE  REVERBERATION  ENVELOPE  arl  • ut 
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'rhe  next  phase  of  the  problem  consisted  of  subjectin^^  the  fluctuating 
component  to  a time  average,  as  indicated  by  Eq.  (ll).  Thus,  if  y^(t) 
is  the  Ith  function  representing  the  fluctuating  component,  then  the 
time  averaged  output  is 


, /-t+T  /2 

Y (t)  = — I ° y (a)dm 


o ■'t-T  /2 
o 


(16) 


However,  since  y(t)  is  represented  by  digital  samnles  uniformly  spaced 
in  time,  the  output  of  the  time  averager  is 


m=n+N  /2 

.o 


r\ 


(IT) 


m=n-N  /2 
o 


The  variable  n corresponds  to  the  time  t in  Eq.  (l6),  whereas  the  number 
of  samples  is  associated  with  the  integration  time  T^.  An  example 
of  the  time  average  is  shown  in  Fig.  3.  The  averaging  time  T^  is  1 msec, 
which  is  the  transmitted  signal  length. 
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TIME  AVERAGE  OF  FLUCTUATING  COMPONENT 


FIGURE  3 

EXAMPLE  OF  THE  TIME  AVERAGE  OF  THE  FLUCTUATING  COMPONENT 
OF  THE  ENVELOPE  FUNCTION  OF  REVERBERATION 
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V.  ESTIMATION  OF  THE  INTENSITY  SPECTRIM 


The  data  analyzed  in  this  study  are  necessarily  of  finite  duration 
and,  therefore,  we  do  not  have  an  infinitely  long,  stationary  random 
prix:ess  available  for  the  estimation  of  the  intensity  spectrum.  Thus, 
we  begin  by  assuming  that  a finite  sized  ensemble  of  functions  representa- 
tive of  the  stochastic  process  (in  this  case  the  fluctuating  component 
of  the  envelope)  is  available  with  each  function  defined  over  the 
interval  (t^,!^).  The  total  time  span  of  our  digital  data  is  the 
interval  (0,  130)  msec,  as  Fig.  1 indicates. 

The  first  step  in  the  estimation  of  the  intensity  spectrum  is  the 
computation  of  the  Fourier  transform  of  each  digital  record.  This  gives 


yAt) 


-iut 


dt 


o)  = 2Ttf  , (18) 


where  thr  index  i vsiries  from  1 to  150,  and  Interval 

chosen  for  the  analysis.  Although  Eq.  (l8)  represents  the  transform  of 
a continuous  function,  the  computation  is  performed  with  digital  samples 
resulting  in  a discrete  Fourier  transform.  The  computation  tool  is  the 
Fast  Fourier  Transform.^ 


The  estimate  of  the  intensity  spectr\jm  is  an  ensemble  average. 


W (f) 

y 


U-f  I 

150  (t  -t.) 
1«1  2 1 


(19) 


where  “ (39,  65)  msec  (see  Fig.  l). 
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The  spectral  estimates  of  both  the  fluctuating  component  and 
its  time  average  (Eq.  (ll))  are  shown  in  Fig.  U.  The  estimates, 
which  are  compared  to  the  appropriate  theoretical  expressions  (Eq.  (7) 
and  Eq.  (ll)),  are  plotted  in  decibels  as  a function  of  frequency  over  the 
interval  (O.O,  2.0)  kHz.  To  allow  a comparison  of  the  shapes  of  these 
curves,  the  maximum  values  of  the  experimental  estimates  were  normalized. 
Over  the  interval  (0.2,  2.0)  kHz  the  experimental  spectra  agree  quite 
well  with  the  theoretically  predicted  values.  Disagreement,  however, 
is  noticeable  in  the  interval  (0.0,  0.2)  kHz.  Here  the  experimental 
spectra  are  approximately  2 or  3 dB  greater  than  predicted.  The  reeison 
for  this  lack  of  symmetry  is  not  presently  known,  but  other  data, 
generated  Independently  of  the  data  analyzed  here,  will  be  studied  in 
a similar  manner  to  determine  if  this  disagreement  still  exists  at  the 
lower  frequencies. 


RELATIVE  LEVEL  OF  INTENSITY  SPECTRA  - dB 


FIGURE  4 

INTENSITY  SPECTRA  OF  FLUCTUATING  COMPONENT 
OF  REVERBERATION  ENVELOPE  (o). 

TIME  AVERAGE  ("SMOOTHED")  OF  FLUCTUATING 
COMPONENT  OF  REVERBERATION  ENVELOPE  (b). 
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VT.  RESEARCH  PLANNED  FOR  THE  NEXT  QUARTER 

Other  statistical  properties  of  the  reverberation  envelope  will 
be  investipated.  For  example,  an  important  assumption  in  the  theoretical 
development  is  made  in  Eq.  (3),  which  relates  the  covariance  of  the 
narrowband  reverberation  to  the  covariance  of  the  envelope.  This  assump- 
tion will  be  tested  with  the  experimental  data,  and  the  results  will 
likely  appear  in  the  next  progress  report. 

Additional  reverberation  data  have  been  generated  at  LaJee  Travis 
>U)d  recorded  on  magnetic  tape.  These  new  data,  which  contain  FM  in 
addition  to  cw  reverberation,  will  be  digitized  and  analyzed  to  determine 
their  covariance  and  spectral  structure.  Hopefully  the  first  results 
of  the  effects  of  time  averaging  deterministic  signals  corrupted  by 
the  additive  reverberation  will  be  reported. 
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